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Abstract 

We study instantons in QCD with many colors. We first discuss a number of qualitative argu- 
ments concerning the large N c scaling behavior of a random instanton ensemble. We show that 
most hadronic observables are compatible with standard large N c counting rules provided the av- 
erage instanton size is 0(1) and the instanton density is 0(N C ) in the large N c limit. This is not 
the case for the topological susceptibility and the mass of the n' ■ For these observables consis- 
tency with conventional large N c counting requires that fluctuations in the instanton liquid are 
suppressed compared to Poissonian fluctuations. Using mean field estimates and numerical simu- 
lations we show that the required scaling behavior of the instanton density is natural in models in 
which the instanton density is regularized in terms of a classical repulsive core. We also show that 
in these models fluctuations of the topological charge are suppressed and that m 2 , = 0(1/N C ). We 
conclude that the instanton liquid model is not necessarily in conflict with the 1/N C expansion. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD), the theory of the strong interactions, is an essentially 
parameter free theory. The QCD lagrangian contains a coupling constant, g, but because of 
the phenomenon of dimensional transmutation this parameter is replaced by a dimensionful 
scale parameter, Aqcd- As a consequence, there is no obvious expansion parameter in 
QCD. 't Hooft proposed to use the number of colors, N c , as a parameter |]J. He suggested 
that QCD simplifies in the limit N c — > oo and that 1/N C can be treated as an expansion 
parameter. 

We do not know how to solve QCD in the large N c limit. However, by analyzing the 
N c dependence of classes of Feynman diagrams we can make certain qualitative statements 
about the structure of large iV c QCD. We expect, for example, that the masses of mesons and 
glueballs are 0(1) whereas the masses of baryons are 0(N C ). Also, meson decay constants 
are 0(iV c 1//2 ) and meson- meson scattering amplitudes are 0(1/N C ). As a result, large N c 
QCD is a theory of weakly interacting mesons and glueballs. 

The fate of the U(1)a anomaly in the large N c limit is an interesting question. The 
problem is complicated by the fact that if the 9 term, 
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is added to the QCD lagrangian then in perturbation theory there is no dependence on the 
parameter 8. Witten suggested that non-perturbative effects generate ^-dependence in the 
pure gauge theory and that the topological susceptibility, 

d 2 E 



Xtop — , (2) 

is 0(1) in the large N c limit [Q, [|. Standard large N c counting suggests that the contribution 
of fermions to Xtop is subleading in 1/N C . We know, however, that this is not correct, because 
there is no ^-dependence in QCD with massless fermions. Witten argued that this apparent 
contradiction can be resolved if the mass of the rj' meson scales as N~ 1 ^ 2 in the large N c 
limit. Witten Q and Veneziano [|5| derived a relation between the mass of the 7/ and the 
topological susceptibility in pure gauge theory, 
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"V = Xtop- (3) 
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Using xtop = 0(1) and f% = 0(N C ) we observe that indeed mi, = 0(1/N C ). This result 
implies that the U(1)a anomaly is effectively restored in the large N c limit. 

The ^-dependence of vacuum energy is related to topological properties of QCD. In the 
semi-classical approximation these features can be described in terms of instantons. Instan- 
tons are localized field configurations that carry topological charge || 



Qt 0P = ^g;g; = ±i. (4) 

In 1978 Witten observed that classical effects, such as instantons, scale as exp(— l/g 2 ) ~ 
exp(— iV c ) which seems to contradict the assumption Xtop — 0(1) ||. On the other hand, 
it was also found that the assumption that topological fluctuations are semi-classical leads 
to a very successful picture of the QCD vacuum, termed the instanton liquid model |7], ||, 
|9], [Hj. The instanton liquid model postulates that for N c = 3 the density of instantons 
is approximately (N/V) ~ 1 fm~ 4 while the average size is p ~ 1/3 fm. These numbers 
reproduce the topological susceptibility in the pure gauge theory Xtop — 

(200 MeV) 4 and the 

chiral condensate (ipip) — —(230 MeV) 3 . More detailed calculations show that the instanton 
liquid model successfully describes an impressive amount of data on hadronic correlation 
functions 0,0. 

Topological properties of the QCD vacuum have also been studied in lattice QCD. It was 
found that the topological susceptibility in pure gauge QCD is Xto P — (200 MeV) 4 P"3"| , as 
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predicted by the Witten- Veneziano relation equ. (£J). However, it was also observed that 
the topological susceptibility is very stable under cooling, and appears to be dominated by 
semi-classical configurations. Lattice simulations also seem to confirm the values of the key 
parameters of the instanton liquid jl4| , (N/V) ~ 1 fm -4 and p ~ 1/3 fm. Recent work has 
focused on the relation of instantons and low-lying eigenstates of the Dirac operator. The 
instanton model predicts that the lowest eigenstates of the Dirac operator, which dominate 
chiral symmetry breaking, are linear combinations of localized, approximately chiral states 



associated with the fermionic zero modes of individual instantons and anti-instantons |T5 
This picture has been confirmed by lattice calculations |l6l IT? . 



In light of these results the question arises whether the success of the instanton model 
can be reconciled with the large N c expansion. In order to address this problem we have 
investigated the predictions of the instanton model for QCD with many colors. This paper 
is organized as follows. In Sect. [TJ we present qualitative arguments based on a random 



instanton liquid. In Sects. [Hj and |V| we discuss analytic results obtained in the mean field 
approximation. In Sects. [VJ and [VJ] we discuss the results of numerical simulations of the 
interacting instanton liquid. 

II. RANDOM INSTANTON LIQUID 

The simplest picture of the instanton liquid in QCD is the random instanton liquid model 
proposed by Shuryak ||. This model is based on the assumption that the instanton ensemble 
is characterized by two key parameters, the instanton density (N/V) and the typical instan- 
ton size p. Except for the size, the collective coordinates of the instantons are distributed 
randomly. The model refers to an instanton liquid rather than an instanton gas because 
interactions between instanton cannot be ignored completely. This is particularly clear if 
massless fermions are present. If interactions are ignored then every instanton contributes 
an exact zero mode to the spectrum of the Dirac operator and the total density of instantons 
is exactly zero. Interactions lead to mixing between the zero modes and qualitatively change 
the spectrum of the Dirac operator. As a consequence, chiral symmetry is spontaneously 
broken and the density of instantons is non-zero. 

In order to understand the large N c behavior of the instanton liquid we have to understand 
the scaling behavior of the two parameters introduced above. In this paper we shall argue 
that 

y = 0(N c ), p = 0{\). (5) 

In section [ITT| we will show that equ. (|5[) follows from the partition function of the instanton 
liquid if the instanton size distribution is regularized by a classical repulsive core. Before we 
discuss this result we would like to provide some plausibility arguments and explore simple 
consequences of the scaling relations. 

We first consider the average instanton size p. At one-loop order the instanton action 
is given by So = (87c 2 )/g 2 = — Mog(pA) where b = (HN c )/3 is the first coefficient of the 
beta function in pure gauge QCD. In the 't Hooft limit N c —>■ oo with g 2 N c = const we 
expect Sq = 0(N C ) and p = 0(1). We have to keep in mind, however, that instantons 
in QCD come in all sizes. In the semi-classical approximation we cannot study instantons 
with size p ~ A -1 and action So ~ 1. There is evidence, however, that chiral symmetry 
breaking in QCD with N c = 3 colors is dominated by small instantons. A priori it is not 
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clear whether this remains true in the large N c limit. Our strategy in this work is to assume 
that the semi-classical approximation remains valid and to show that this assumption leads 
to a consistent picture. 

In the random instanton liquid model the instanton density is related to the non- 
perturbative gluon condensate 

Standard N c counting suggests that (g 2 G 2 ) = 0(N C ) and we are lead to the conclusion that 
(N/V) = 0(N C ). This is also consistent with the expected scaling of the vacuum energy. 
Using equ. (§) and the trace anomaly relation 

b 

the vacuum energy density is given by 

b /N" 



( T w) = -^(9 2 G%G°), (7) 



(8) 



Using (N/V) = 0(N C ) we find that the vacuum energy scales as e = 0(N 2 ) which agrees 
with our expectations for a system with N 2 gluonic degrees of freedom. We should note that 
equ. (H) requires interactions between instantons. The partition function of a completely 
non-interacting gas of instantons gives e ~ (N/V) ~ N c . 

We also note that (N/V) = 0(N C ) implies that the effective "packing fraction" of instan- 
tons remains constant in the large N c limit. Instantons with topological charge Qt op = ±1 in 
SU(N C ) QCD are embeddings of SU(2) instantons. Since the number of mutually commut- 
ing SU(2) subgroups of SU(N C ) scales as N c we can have 0(N C ) instantons which overlap 
in space but are in fact weakly interacting. Witten argued that the only alternatives are 
S = 0(N C ) and (N/V) = 0(e~ Nc ) or S = 0(1) and (N/V) = 0(e Nc ). However, as we 
shall see in the next section, it is possible for the instanton density to remain finite even if 
Sq = 0(N C ) because the large entropy of instantons in SU(N C ) can overcome the exponential 
suppression due to the action. 

If instantons are distributed randomly then fluctuations in the number of instantons and 
anti-instantons are expected to be Poissonian. This leads to the predictions 

(N 2 )-(N) 2 = (N), (9) 
(Q 2 ) = (N), (10) 



where N = Nj + Na is the total number of instantons and Q = Nj — Na is the topological 
charge. Equ. (110) implies that 

(Q 2 ) (N 



XtO P = — ={y). (11) 

Using (N/V) = 0(N C ) we observe that Xtop — 0(N C ) which is in contradiction to Witten's 
assumption Xtop — 0(1). However, as we shall see in the next section, interactions between 
instantons cannot be ignored in the large N c limit and the fluctuations are suppressed 
compared to equs. (PJlOl). 

Finally, we would like to study chiral symmetry breaking in a random instanton liquid. 
For definiteness, we will consider the case Nf = 2 but the conclusions are of course inde- 
pendent of the number of flavors. Instantons induce an effective 2iVy-fermion lagrangian. 
After averaging over the color orientation of the instanton the effective lagrangian is given 
by MM 



C = J n(p)dp ^^-jr € flf2 e gig2 ( ^ N ^ N 1 ^LJ^R, 9 ,){^LJ 2 ^R, g2 ) (12) 

We observe that the explicit N c dependence is given by l/N%- This is again related to the 
fact that instantons are SU(2) objects. Quarks can only interact via instanton zero modes 
if they overlap with the color wave function of the instanton. As a result, the probability 
that two quarks with arbitrary color propagating in the background field of an instanton 
interact is 0(l/iV 2 ). 

Chiral symmetry breaking can be studied in the mean field approximation. We will 
address this problem in much more detail in the following section but we can give a simple 
qualitative argument here. In the mean field approximation we can derive a gap equation 
for the spontaneously generated constituent quark mass. The gap equation is of the form 

M = GN C ( d \ 9 M , 9 , (13) 
J (2tt) 4 M 2 + P' v ; 



where M is the constituent mass and G is the effective coupling constant in equ. (12). The 



factor iV c comes from doing the trace over the quark propagator. The coupling constant 
G scales as 1/N C because the density of instantons is 0(N C ) and the effective lagrangian 
contains an explicit factor 1/A^ 2 . We conclude that the coefficient in the gap equation is 
0(1) and that the dynamically generated quark mass is 0(1) also. This also implies that 
the quark condensate, which involves an extra sum over color, is 0(N C ). 



III. THE MEAN FIELD APPROXIMATION 



In this section we shall make the arguments presented in the previous section more 
quantitative. Instead of considering the average instanton size and instanton size density to 
be arbitrary parameters, we shall determine these quantities from the partition function of 
the instanton ensemble. This means, in particular, that both p and (N/V) are expressed in 
terms of the fundamental scale parameter of QCD. For this purpose it is essential to take 
the interaction between instantons into account. If instantons are semi-classical, Si ns t ^> 1, 
and if the interaction between instantons is weak, Si nt -C Si nst , this can be accomplished 
using mean field methods || ED|, ^T|]. In this section, we shall follow the variational method 
of Diakonov and Petrov 0. 

We consider the partition function for a system of instantons in pure gauge theory 

, v 7 +v A 

Z = Jm^\ II JWin{ Pl )}e W {~S mt ). (14) 
Here, Qj = (zj,pj,Uj) are the collective coordinates of the instanton / and n(p) is the 



semi-classical instanton distribution function 22 



n(p) = C Nc I — I p exp 
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9(p) 



(15) 



0.466 exp(-1.679iV c ) 
Cn ° ~ (N e -l)\(N e -2)\ ' (16) 

-Mog(pA), b=^N c . (17) 
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We have denoted the classical instanton interaction by S% n t. If the instanton ensemble is 
sufficiently dilute we can approximate the instanton interaction as a sum of two-body terms, 
Si n t = J2u Si j- F° r a we U separated instanton- anti- instanton pair the interaction has the 
dipole structure 

^ = -^!#M 2 (l-4cos 2 #). (18) 
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Here pi t A are instanton radii and Rja is the instanton-anti-instanton separation. The relative 
color orientation is characterized by a complex four- vector = 7^tr(L r / J 4T+), where Via — 
UiU A depends on the rigid gauge transformations that describe the color orientation of the 
individual instanton and anti- instanton and r+ = (r, —i). We have also defined the relative 
color angle cos 2 6 = \u ■ R\ 2 /\u\ 2 . The dipole interaction is valid if Rj A ^> Pi Pa- We will 
specify the interaction at shorter distances below. 



Diakonov and Petrov suggested to analyze the partition function equ. (|14]) using a vari- 
ational single-instanton distribution p(p) ||. The corresponding partition function is 

i Nj+Na r i 

^wai\ n / d!i '"M = W^°> N ' +A " (19) 

where po = J dpp(p). The exact partition function is 

Z = Z 1 (exp(-(S-5 1 ))), (20) 

where 5* is the full action, Si = log(/x(p)) is the variational estimate and the average (.) is 
computed using the variational distribution function. The partition function satisfies the 
bound 

Z>Ziexp(-(S-Si)), (21) 

which follows from convexity. The optimal distribution function p(p) is determined from 
a variational principle, (5 log Z)/ (5 p(p)) = 0, where Z is computed from equ. (pl|). One 
can show that the variational result for the free energy F = — \og(Z)/V provides an upper 
bound on the true free energy. 

The calculation of (S — Si) reduces to the calculation of the average instanton interaction 
{Si n t). Since the variational ansatz does not include any correlations, we only need to average 
the instanton interaction over the collective coordinates of the two instantons. The dipole 
interaction (|IBD vanishes when averaged over all color orientations. In || it was proposed to 
compute the instanton interaction at all distances using a specific ansatz (called the "sum 
ansatz" ) for the two-instanton configuration. The result is that both the instanton-instanton 
(//) and instanton-anti-instanton (I A) are repulsive on the average. We find 

I O \ ^ 2 2 2 2 27 N c 2 /oo\ 

(Sto) = -Jl PjPj, 7 = X^vf^T 7 " ' (22) 

The interaction contains an explicit factor N c / (N^ — 1) ~ 1/N C which reflects the probability 
that two random instantons overlap in color space. Since the classical action scales as 
So ~ l/g 2 we find that the average interaction between any two instantons is 0(1). Applying 
the variational principle, one finds 1| 



Mp) = ™(p)exp 



(23) 



V 

where (3 = (3(j>) is the average instanton action and p 2 is the average size. We observe that 
the single instanton distribution is cut off at large sizes by the average instanton repulsion. 
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The instanton density and average size are given by 



7"* 



c N j 2N ^{v){^r^y l \ (24) 



vV \ 1 ' 2 6-4 



" 2 = {-m • v= — (25) 



(26) 



We note that A is the only dimensionful parameter. The free energy is given by 

which is in agreement with the trace anomaly. We can now study the dependence of (N/V) 
and p on N c , see Fig. |I[ We note that to one- loop order the scale in the pre-exponent (3(p) 
is not well determined. In practice we assume that j3 = N c so with s = 5. Changing so does 
affect both (N/V) and p but the main effect can be absorbed in the scale parameter. The 
remaining dependence on sq is very weak. 

Fig. [I] shows that for N c > 4 the average instanton size is essentially constant while 
the instanton density grows linearly with N c . This is easily verified by inspecting equ. (p4|) . 
Expanding log(N/V) in powers of N c and log(iV c ) we observe that independent of the details 
of the interaction the instanton density scales at most as a power, not an exponential, in 
N c |23j. Using the fact that 7 2 = 0(1/N C ), which is equivalent to (S int ) = 0(1), we find 
that (N/V) = 0(N C ). This result depends on the instanton interaction, but as we noted 
above, 7 2 = 0(1/N C ) is a consequence of the fact that instantons are SU(2) gauge field 
configurations. 

There is a simple argument that explains why the instanton density scales as the number 
of colors. In our model, the size distribution is regularized by the interaction between 
instantons. This means that there has to be a balance between the average single instanton 
action and the average interaction between instantons. If the average instanton action 
satisfies S = 0(N C ) we expect that (Sffi) = 0(N C ) also. Using (S^ t ) = (N/V)(S int ) and 
the fact that the average interaction between any two instantons satisfies (Si nt ) = 0(1) we 
expect that the density grows as iV c . 

Fig. |2| shows the instanton size distribution for different numbers of colors. We observe 
that the number of small instantons is strongly suppressed as N c — > oo but the average 
size stabilizes at a finite value p < A -1 . We also note that there is critical size p* for 
which the number of instantons does not change as iV c — > oo. The value of p* is easy to 
determine analytically. We write n(p) = exp(N c F(p)) with F(p) = alog(p) + bp 2 + c where 



the coefficients a, b, c are independent of A" c in the large N c limit. The critical value of p is 
given by the zero of F(p). We find p* = 0.49A -1 . The existence of a critical instanton size 
for which n(p*) is independent of N c was discussed by [[24], ^5], |26|]. The problem was studied 



on the lattice by Lucini and Teper [B7[, who find p* = 6a = 0.43 fm. 

Next we wish to study fluctuations in the instanton liquid. Fluctuations in the net 
instanton number are related to the second derivative of the free energy with respect to N. 
We find 

(N 2 ) - (iV> 2 = 7 (AO. (27) 
o 



This result is in agreement with a low energy theorem based on broken scale invariance |28 
/ d±x {(s 2 GWG 2 (*)> - (s 2 G 2 (0)) 2 } = l^(9 2 G 2 ). (28) 



(327T 2 ) 2 

This result is very general and based solely on the renormalization group equations. The 
left hand side is given by an integral over the field strength correlator, suitably regularized 
and with the constant term (G 2 ) 2 subtracted. For a dilute system of instantons equ. ( 28|) 



reduces to equ. ([271). The result fl27|) shows that fluctuations of the instanton ensemble are 
suppressed by 1/N C . This is agreement with general arguments showing that fluctuations 
are suppressed in the large A" c limit. We also note that the result ( p7|) clearly shows that 
even if instantons are semi-classical, interactions between instantons are crucial in the large 
A^ limit. 

Fluctuation in the topological charge can be studied by adding a #-term to the partition 
function flTj|). We find 

(Q 2 ) = (AO, (29) 
which is identical to the result in the random instanton liquid and not in agreement with 



Witten's hypothesis Xtop = 0(1). However, Diakonov et al. noticed that equ. (|29f) is a 
consequence of the fact that in the sum ansatz the average interaction between instantons 
of the same charge is identical to the average interaction between instantons of opposite 



charge [29]. In general there is no reason for this to be the case and more sophisticated 



instanton interactions do not have this feature JSD, 51], |32[. If r denotes the ratio of the 
average interaction between instantons of opposite charge and instanton of the same charge, 
r = (S IA )/(S„}, then g| 



« 2 > - V^T) W - (30) 
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This result shows that for any value of r ^ 1 fluctuations in the topological charge are 
suppressed as N c — > oo. We also note that xto P = 0(1), in agreement with Witten's 
hypothesis. 

IV. CHIRAL SYMMETRY BREAKING 

In this section we wish to study chiral symmetry breaking in the mean field approxima- 
tion. This can be done by studying the Dyson-Schwinger equation for the quark propagator 
or by analyzing the spectrum of the Dirac operator. In this section we wish to use the more 
microscopic approach and analyze the spectrum of the Dirac operator. In a basis spanned 
by the individual zero modes of the instantons and anti-instantons the Dirac operator has 
the structure 

m=[ ° T,A ), (3D 

where Tia is the overlap matrix element of the Dirac operators between an instanton and 
anti-instanton zero mode. The matrix elements depend on the collective coordinates of the 
instanton. If the interaction between instantons is weak, the matrix elements are distributed 
randomly with zero average, but the second moment of Tia is non-zero. Averaging over the 
positions and orientations of the instantons we get 



The factor l/N c comes from the average over SU(N C ). Equ. fl32|) implies that the average 
matrix element of the Dirac operators decreases as 1/N C but the second moment in the zero 
modes zone is 0(1). If the matrix elements are distributed according to a Gaussian unitary 
ensemble, the spectral density is a semi-circle 

1/2 

iraV V" 4a 2 



with a 2 = \Tf A \. We observe that the width of the zero mode zone is related to o~, which 
is 0(1). According to the Banks- Casher formula the quark condensate is related to the 
spectral density at zero virtuality 

. , 1 /3A C A\ 1/2 

{ ^--- p {— v) ■ (34) 
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Because (N/V) = 0(N C ) the quark condensate also grows as N c . This, of course, is the 
expected behavior. We note, however, that the linear growth in N c is really a combination 
of two effects. The linear increase in the number of modes N provides one factor of N^ 2 
and the decrease in the average matrix element \Tja\ contributes another factor of N^ 2 . 

The true spectral density of the Dirac operator in not given by a semi-circle. Schemati- 
cally, the spectrum of the Dirac operator in quenched QCD is shown in Fig. ^. There are 
several notable features. First, in a finite volume there is a certain number of exact zero 
modes. This number is proportional to (xtopV) 1 ^ 2 and therefore scales as N®. We note that 
in the infinite volume limit exact zero modes are not important because their number scales 
as \/V while the total number of states increases linearly with the volume. The second 
feature of the spectrum is the logarithmic enhancement of the spectrum at small virtual- 
ly. This enhancement is an artifact of the quenched approximation. From quenched chiral 
perturbation theory we expect |33 



Here, £ is a parameter that corresponds to the (negative) unquenched chiral condensate 
and mo is the mass of the quenched ghost pole. This mass corresponds to the mass of the 
7/ in the chiral limit of full QCD. Using m\ = 0(1/N C ) we find that the coefficient of the 
logarithmic enhancement vanishes as 1/N C in the large N c limit. This, of course, is consistent 
with the idea that the fermion determinant is not important in the large N c limit. The third 
component of the spectrum is given by the almost zero modes related to chiral symmetry 
breaking. This part of the spectrum is expected to scale as N c , as is the bulk of the spectrum 
which is not related to chiral symmetry breaking. 

At finite N c the spectrum of the Dirac operator in full QCD with light fermions has a 
different behavior. The number of exact zero modes is again proportional to (xtopV) 1 ^ 2 but 
in full QCD the topological susceptibility is suppressed, Xtop — —m{ipip). This implies that 
as N c increases the number of exact zero modes initially increases as N^ 2 and then saturates 
at the value corresponding to the quenched topological susceptibility. In the infinite volume 
limit the spectrum near the origin is linear. Chiral perturbation theory predicts the slope 
of the spectrum pi| 

£ ( (iV? — 4)E, , 1 
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We observe that the constant part of the spectrum grows as N c whereas the linear part is 
independent of N c . This implies that in the large N c limit the spectrum at the origin is flat 
for any number of flavors as long as Nf is not of the order iV c . 



V. THE INTERACTING INSTANTON LIQUID 

In this section we wish to go beyond the mean field approximation and study the par- 
tition function of the instanton liquid using numerical simulations. These simulations take 
into account all correlations between instantons. The numerical techniques are described 



in detail in [35, 36|. In order to perform these simulations we have to fully specify the 



instanton interaction. We have used the "streamline" interaction determined in |3T], |32 



The streamline solution is characterized by the fact that the action of the approximate 
instanton-anti-instanton solution is a local minimum except in the direction of the "valley" 
in configuration space that connects a well separated IA pair with a very close pair. There 
is no interaction between two instantons of the same charge. The interaction between two 
instantons of opposite charge approaches — 25o if the relative color orientation is attractive. 

This implies that the streamline interaction lacks the repulsive core that is required to 
stabilize the instanton ensemble at the classical level. In order to correct this problem we 
have added a purely phenomenological core to the instanton interaction. The interaction is 
given by 

o 8 ^ U 2 ,_ R 2 + p 2 i + p 2 a | / (R 2 + p] + p 2 a) 2 X /2 ^ 

5core -^ N ' A - 2 PlPA + { 4M l ) (3?) 

in both II and I A channels. The dimensionless parameter A controls the strength of the 
core and is adjusted to reproduce the phenomenological diluteness p 4 (N/V) of the instanton 
ensemble. In our simulations we have used A = 128. We note that the hard core interaction 
equ. (|3"7j) only acts between instantons that overlap in color space. As a result, we expect 
(Score) — 0(1), in agreement with the interaction used in the mean field treatment. 

One might argue that there should not be any interaction between instantons of the same 
charge because there is a 2 x (2A C ) parameter family of two-instanton solutions with action 
2 x (8n 2 )/g 2 [38]. However, there are two phenomena that lead to an effective instanton- 



instanton interaction. The first is related to the fact that the collective coordinate measure 
for two close instantons is not just the product of two single instantons measures. Carter and 
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Shuryak argued that this will lead to an effective short range instanton-instanton repulsion 
|39|, [4(|. Whether this effect has the same N c dependence as the classical interaction equ. fl3"7|) 
is an important problem. The second effect is that quantum corrections to a charge two 
instanton solution do not factorize. As a quantum correction, this effect is naively suppressed 
by 1/Sq ~ 1/iVc, but the suppression can be overcome if instantons with different color 
orientation interact. 

In order to determine the instanton density and the free energy of the instanton liquid we 
have to compute the partition function equ. (|14|) . Monte Carlo simulations are ideally suited 
for computing expectation values, but they do not directly provide the partition function. 
It is possible, however, to compute free energy differences. This implies that Monte Carlo 
simulations can be used to compute the ratio of two partition functions. In practice we 
calculate the ratio of the exact partition function and the variational estimate equ. (19). 
The exact partition function is given by |36| 



\ogZ = \og{Z 1 )- C da{{S-S 1 )) a . (38) 
Jo 

Here, Z\ is the variational partition function equ. ( |T9| ) and S\ is the variational action. The 
expectation value (.)„ is determined using the interpolating action S a — S\ + a(S — Si). S a 
reduces to the variational action for a = and the exact action for a = 1. 

The free energy energy of the instanton liquid as a function of the instanton density for 
N c — 3, . . . , 6 is shown in Fig. |4j. The equilibrium density is determined by the minimum 
of the function F(N/V). The dependence of equilibrium density and the free energy on 
the number of colors is shown in Fig. || We observe that (N/V) increases linearly with N c 
whereas the free energy is quadratic. The slope of (N/V) as a function of N c is small, in 
agreement with the mean field result shown in Fig. |l[ In contrast to the mean field result 
the linear behavior already sets in at small N c ~ 3. We have checked the stability of our 
results to changing the strength of the hard core interaction and including higher order 
corrections in the QCD beta function. Both changes affect the results quantitatively but 
not qualitatively. However, the results are crucially dependent on the assumption that the 
parameter A in equ. ([37]) is not a function of N c . 

We have also studied the instanton size distribution, the topological susceptibility and 
the spectrum of the Dirac operator for different numbers of colors. In order to be able to 
distinguish more clearly between different scenarios we have not used the exact instanton 
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density determined in Fig. |] but have simply scaled (N/V) ~ iV c . At large N c this will 
only introduce errors that are suppressed by 1/N C . The instanton size distribution is shown 
in Fig. H As expected small instantons are strongly suppressed as the number of colors 
increases. We observe a clear fixed point in the size distribution at p*A ~ 0.27. 

Our simulations were carried out in the total topological charge Q top = sector of the 
theory. We can nevertheless determine the topological susceptibility by measuring the av- 
erage Ql in a sub- volume V 3 x Z 4 of the euclidean box V 3 x L 4 f35|. The finite volume 
susceptibility is given by 

= ^ff(i -£)"'• < 39 > 

The factor (1 — U/L^) -1 takes into account the constraint from overall charge neutrality. 
This correction factor is derived under the assumption that the fluctuations are Gaussian. 
In an ideal calculation L 4 3> U and the correction for overall neutrality is small. The topo- 
logical susceptibilities are shown in Fig. [7[ We observe that Xto P (U) tends to a constant as 
1 4 increases. We identify this constant with the susceptibility in the thermodynamic limit. 
We find that for N c = 3 the topological susceptibility agrees well with the expectation based 
on Poissonian statistics, Xto P — {N/V). For N c > 3, however, fluctuations are significantly 
suppressed and the topological susceptibility increases more slowly than the density of in- 
stantons. Fig. [| shows that our results are consistent with a scenario in which Xtop remains 
finite as iV c -> oo. 

In Fig. ^ we show the spectrum of the Dirac operator for N c = 3, . . . , 6. Since we calculate 
in the Q top = sector of the theory there are no exact zero modes. We clearly observe the 
enhancement of the spectral density near A = 0, but we also note that this enhancement 
becomes weaker as the number of colors increases. The chiral condensate for m q = 0.1 A is 
shown in Fig. |8|. We clearly see that (qq) is linear in N c . 

In the recent literature a number of authors have studied the role of instantons in chiral 
symmetry breaking by analyzing the local chirality X(x) of low lying eigenstates of the Dirac 
operator [T7]]. The quantity X(x) is defined by 

^fa+wHw^T- (40) 

where ip is an eigenfunction of the Dirac operator iJ^ip = \ip in a given gauge configuration. 
In order to study chiral symmetry breaking one only considers the lowest few eigenvec- 
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tors. The instanton liquid model predicts that these eigenvectors are linear combinations 
of instanton and anti-instanton zero modes. This implies that at points x where the wave 
function ijyij) is large it is either left or right handed, X(x) ~ ±1. In order to test this 
prediction one has to choose a cutoff on the eigenvalue A and some cutoff on the magnitude 
of ijj^ip. Typically, the points x are restricted to the top few percent of the eigenfunction. 
The instanton model suggests that this fraction should be no larger than the diluteness of 
the instanton liquid, and that the maximum eigenvalue should be smaller than the width 
\Tja\ of the zero mode zone. 



In Fig. [10] we show numerical results for the local chirality distribution in the instanton 
liquid model. We have included all states in the zero mode zone and used the top 30% of the 
wavefunction. The N c dependence is similar if more restrictive cuts are used. We observe 
that the double peak structure is very pronounced for all N c — 3, . . . , 6. There is a small 
shift of the peaks toward smaller values of X as the number of colors increases. We have 
verified that for N c > 10 the double peak structure disappears. This is related to the fact 
that the instanton density increases and instantons overlap in space. The instanton liquid 
remains dilute, however, because instantons do not overlap in color space. This is easily 
verified by computing chirality distributions for eigenstates of the Dirac operator projected 
on a SU(2) subgroup. Cundy et al. computed local chirality distributions in lattice QCD 
for N c = 2, . . . , 5 In their calculations the double peak structure disappears much more 



quickly as the number of colors increases. This implies that either the instanton ensemble is 
not as dilute as we have assumed, or that mixing with non-zero modes is more important. 
Given the fact the value of p* in our simulations is also smaller than the value obtained by 
Lucini and Teper, it seems likely that the instanton liquid in QCD is not quite as dilute as 
suggested in [|J. 

VI. HADRONIC CORRELATION FUNCTIONS 

We have also studied hadronic correlation functions in the instanton liquid for different 
numbers of colors. We have considered, in particular, correlation functions of currents with 
the quantum numbers of the pion, the rho meson, and the rf meson. The currents are given 
by 

j'tt = di^u, f p = d^u, j v > = -j=(ui-/ 5 u + di^d). (41) 
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The pion and rho meson are flavor non-singlet mesons. The corresponding correlation func- 
tions only involve the connected contribution 



n /=1 (x) = (Tr s ab (o,x)r\s ba (x,o)r ) 



(42) 



with T = Z75,7 M for the pion and rho meson, respectively. Here, S ab denotes the quark 
propagator in a given gauge configuration, a, b are color indices, and the trace runs over Dirac 
indices. The average (.) is performed with respect to the partition function equ. (|14T) . The rf 
meson is a flavor singlet meson and the correlation function has an additional, disconnected, 
term. We find 



n /=0 (a;) = (Tr S ab (0,x)TS ba (x,0)T ) -2(Tr[S aa (0,0)r] Tr S bb (x,x)T ) 



iba i 



ibb, 



(43) 



with T = 275. The masses of the tt, p and 77' mesons can be extracted by fitting the correlation 
functions to a spectral representation. At large euclidean separation the correlation function 
is determined by the lowest hadronic resonance. A good model for the contribution of excited 
states is given by the free quark- ant i- quark continuum above some threshold invariant mass 
y/so. The threshold roughly corresponds to the location of the first excited state. For the 
pion and rho meson we use 



U^(x) = XlD(m n , x) + — -| / dssD(y/s,x 

07T J so 



SO 

2N r < - 



n p0) = fp m pD{m p ,x) + ^| / dssD(^/s,x), 



(44) 
(45) 



where D(m, x) is the euclidean space propagator for a free scalar meson with mass m 



m 

D{m,x) = —7- 



Ki(mx). 



(46) 



In full QCD the i]' meson can be described by the same spectral representation as the 
pion. In quenched QCD, however, the spectral function in the rj' channel is not positive 
definite. There is a well established method for dealing with this problem. In quenched 
chiral perturbation theory it is assumed that the disconnected part of the 77' correlation 



function corresponds to a ghost pole in the spectrum [|4*2"| , f43| , HJ. In momentum space the 
pole contributions are given by 

A! XI o 1 



IV (<z) 



-VTLr 



q 2 + ml q 2 + ml °q 2 



m 



2 ' 



(47) 
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Going from quenched to unquenched QCD corresponds to summing the geometric series 

M») - f + (48) 

and mi, = m\ + m 2 in full QCD. This result is expected to be exact in the large N c limit. 
In addition to that there is evidence from the lattice that the mass of the rj' meson can be 



extracted from the ghost pole propagator even for N c = 3 [33J]. In coordinate space equ. ( 4"?D 
corresponds to 



with 



U^(x) = \lD(m in x) - X^mlD'irriT^x) + — / dssD(y/s,x), (49) 

07T Z 7 Sn 



D'( m ,i) = ^oH. (50) 



We have used the spectral representation equ. to determine the quenched rf mass mo 
for different N c . 

Before we come to the numerical studies we would like to present several analytical results. 
If we study correlation functions at distances small compared to the average separation 
between instantons, x <C (iV/V) -1 / 4 , it is sufficient to take into account the contribution 
from the closest instanton only. The quark propagator is given by 

S( X ,y)* Mx -'M<y-') (51) 
m* 

where ipo(x) is the zero mode wave function and z is the location of the instanton. The effect 
of all other instantons only enters via the effective mass m* ~ np(2/N c ) l / 2 (N /V) 1 ' 2 . The 
single instanton contribution to the correlation function is found by inserting the propagator 



equ. (pID into the equations ( ^4) and ( |43|) and averaging over the collective coordinates of 
the instanton. 

Because of the chiral structure of the rho meson current there is no zero mode contribution 



to the rho meson correlation function. In the pion and 77' meson channel we find [[45 

n- W = ± /^ p) ^J^ j|f (_^_ + f logi±|)}, ,52, 

with £ 2 = x 2 /(x 2 + 4p 2 ). As is well known, the instanton contribution is attractive in the 
pion channel, and repulsive in the 77' channel. The contribution of one instanton of given 
size does not involve any factors of iV c , and is the same, up to the sign, in the pion and 
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rf channel. This is illustrated in Fig. [□]. Perturbative contributions to the disconnected 
correlation function are suppressed by a factor 1/N C compared to the connected correlator. 
However, the single instanton contribution to the disconnected correlator is exactly the 
same as the instanton contribution to the connected correlator. After integration over the 
instanton distribution the correlation function scales as N c because the instanton density is 
proportional to N c . This is the expected behavior in the case of the pion correlation function 
but it implies that in the single instanton approximation the anomalous tt — rf splitting does 
not disappear in the large N c limit. 

In order to go to large distances, x > (N/V)^ 1 ^ we have to resum the instanton inter- 



action. This can be achieved using the mean field (Hartree) and RPA methods [46|. The 
mean field approximation gives the constituent quark propagator 

Sn(x) - [*E-<#* t> + iM{p) 

where M(p) is the dynamically generated quark mass. The momentum dependence of M is 
determined by the Fourier transform of zero mode wave function. In the rho meson channel 
there is no direct instanton induced interaction and the correlation function is given by two 
non-interacting constituent quarks. We have 

U r (x) MFA = N c Tr [S Q (x)TS Q (-x)T] , (54) 

where the trace is over the Dirac indices only and T = 7 M . In the pion and rf channel there is 
a direct instanton interaction that can be resummed using the RPA. We have fTO, pi, H7L H8 



u^(x) = nr A (x) + n^/ (55) 

with 

KA X ) = N < (^f) J *qe** r 5 (g) T ^^yr 5 (g). (56) 
The loop and vertex functions C5 and T 5 are given by 

CM = 4iV c y J — (Mi2+pf)(M|+pi) , (57) 
rM r d'p {M x M 2 fl\M x M 2 - Vl -p 2 ) 

i5W " J (2tt)* (M?+pl)(Mi+pl) ' [ ™ } 
where p\ = p + q/2, p 2 = p — q/2 and Mi j2 = M(p lj2 ). Using (N/V) ~ N c and M ~ 1 we 
observe that the pion and rf correlation functions in the large -/V c limit depend on N c only 
through an overall factor of N c . This implies f% = 0(N C ) and m w , = 0(1). 
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Numerical results for the constituent quark mass Mq as well as the pion mass m w and 
pion decay constant f n are shown in Fig. [12]. We have used the average instanton size and 
density determined in section |T|. The current quark mass is m q = 0.025A. There is some 
variation in the constituent quark and pion masses for small iV c < 10 but the size of 1/N C 
corrections is not large, about 20% for N c = 3. The pion decay constant shows the expected 
N^ 2 behavior but in this case 1/N C corrections are large, about 80% for N c = 3. 

In Fig. ^ we show the correlation functions in the pion, rho meson, and rf meson channel. 
The correlation functions are normalized to free field behavior. The overall factor of N c drops 
out if the correlator is normalized in this way. We observe that the p meson correlation 
function is essentially independent of N c already for small N c . There are substantial 1/N C 
corrections in the 7r and rf channel. The splitting between the 7r and rf correlation functions 
is reduced in going from N c — 3 to N c — 6 but it remains finite and large as N c — ► oo. 

The correlation functions measured in numerical simulations of the instanton liquid for 
N c = 3, ... 6 are shown in Fig. 14. The meson masses extracted from the spectral repre- 
sentation equs. (|44]^5| , |49|) are shown in Fig. |H| The results were obtained from simulations 
with iV = 128 instantons in a euclidean volume VA 4 = V3 x 5.76. V3 was adjusted such that 
(N/V) = (N C /3)A 4 . In order to avoid finite volume artifacts the current quark mass was 
taken to be rather large, m q = 0.2A. We observe that the rho meson correlation function 
exhibits almost perfect scaling with N c and as a result the rho meson mass is practically 
independent of N c . The scaling is not as good in the case of the pion. As a consequence 
there is some variation in the pion mass. However, as one can see from the fit shown in 



Fig. [15|, this effect is consistent with 1/N C corrections that amount to about 40% of the pion 
mass for N c = 3. Finally, we study the behavior of the rf correlation function. There is a 
clear tendency toward U(1)a restoration, but the correlation function is still very repulsive 
for N c = 6. As one can see from Fig. [15| the result is consistent with nit, ~ 1/N C although 
the error bars are quite large. 

For comparison, we show the expected behavior of the correlation functions based on 



standard large N c counting in Fig. |16[ We have used the spectral representation equs. ([4^f49D 
together with the phenomenological values m n = 139 MeV, X n = (450 MeV) 2 , E = 1.3 GeV 
and m = 900 MeV (N c = 3). We assume that A 2 ~ and m\ ~ 1/N C . We observe that 
the rf correlation function only approaches the pion correlation for fairly large values of A^ c . 
For example, the rj' correlation function does not show intermediate range attraction unless 
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N c > 15. The variation in going from N c = 3 to iV c = 6 is not dramatic, in agreement with 
the results shown in Fig. ITJ. 



VII. SUMMARY 

In summary we have studied instantons in the large N c limit of QCD. We have argued 
that it is possible for the instanton liquid model to have a smooth large N c limit which is in 
agreement with scaling relations derived from Feynman diagrams. In this limit the density 
of instantons grows as N c whereas the typical instanton size remains finite. Interactions 
between instanton are important and suppress fluctuations of the topological charge. As 
a result the U(1)a anomaly is effectively restored even though the number of instantons 
increases. Using mean field arguments [[| and numerical simulations we have shown that 
this scenario does not require fine tuning. It arises naturally if the instanton ensemble is 
stabilized by a classical repulsive core. In this case we obtain a picture in which the instanton 
density is large but the instanton liquid remains dilute because instantons are not strongly 
overlapping in color space. Further investigations will have to show whether this scenario 
is indeed correct. For example, it would be useful to study the exact moduli space for 
multi-instanton configurations in the large N c limit. 

Of course, there is no a priori reason why instantons have to be compatible with standard 
large N c counting rules. Instantons are not part of the diagrammatic expansion and do not 
need to satisfy scaling relations derived from diagrams. On the other hand it would be 
hard to reconcile the quantitative success of the instanton liquid model in describing chiral 
symmetry breaking and the mass of the r( with the phenomenological success of the 1/N C 
expansion if instanton effects strongly violate N c counting rules. 

We should emphasize that the numerical results presented in this work only cover fairly 
small values of N c , N c < 6, and that all the analytical results were obtained in the mean 
field approximation. We cannot exclude the possibility that there is a phase transition as 
the number of colors becomes large E|, EM ESI . Carter and Shuryak suggested, for example, 



that for clusters involving 0(N C ) instantons 1/N C suppressed color singlet exchanges become 
dominant and lead to the formation of tightly bound molecules. We did not observe this 
phenomenon in our simulations even if a short range color singlet interaction was included, 
but the number of colors (N c < 10) may have been too small. We also did not investigate 
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the possibility that instantons in the large N c limit melt or dissociate into constituents with 
fractional topological charge. The latter scenario was investigated in the case of the CP^ -1 
model in . 

Our results can be compared to the lattice results of Lucini and Teper |27| and Cundy, 
Teper and Wenger Lucini and Teper find a fixed point in the instanton size distribution, 
in agreement with our results shown in Figs. |2] and |5|. However, they do not find any 
suppression of large size instantons. Cundy et al. studied the local chirality distribution. 
They find that the double peak structure of this distribution disappears in the large N c limit. 
Our results shown in Fig. [lOj are in qualitative, but not in quantitative agreement with their 
findings. In our calculations the effect is much less pronounced. This implies that either the 
instanton liquid is not as dilute as it is in our calculations, that mixing with non-zero modes 
becomes more important as N c increases, or that instantons are no longer semi-classical. 
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FIG. 1: Average instanton size p, instanton density (N/V) and quark condensate (qq) for differ- 
ent numbers of colors N c . The results shown in this figure were obtained using the mean field 
approximation. All quantities are given in units of the QCD scale parameter A. 
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FIG. 2: Instanton size distribution n(p) for different numbers of colors N c = 3, . . . , 10. The results 
show in this figure were obtained using the mean field approximation. 
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FIG. 3: Schematic behavior of the spectrum of the Dirac operator in quenched QCD with many 
colors. 
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(N/V) [A 4 ] 

FIG. 4: Free energy F of the quenched instanton liquid as a function of the instanton density 
(N/V) for N c = 3, . . . , 6 colors. Both (N/V) and F are given in units of the QCD scale parameter. 
The results shown in this figure were obtained using numerical simulations with N = 32 instantons. 
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c 

FIG. 5: Instanton density (N/V) and free energy F in a pure gauge instanton ensemble for N c = 
3, . . . , 6 colors. Both (N/V) and F are given in units of A 4 where A is the QCD scale parameter. 
The dashed lines show fits of the form a\N^ + a.2N c + 03 (for the free energy F) and a,2N c + 03 (for 
the instanton density N/V). 
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FIG. 6: Instanton size distribution in a pure gauge instanton ensemble for different numbers of 
colors. The results were obtained using numerical simulations with N = 128 instantons. 
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FIG. 7: Finite volume topological susceptibility Xtapi}^) m a pure gauge instanton ensemble for 
different numbers of colors. The results were obtained using numerical simulations with iV = 128 
instantons. 
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FIG. 8: Dependence of the chiral condensate (ipip) and the topological susceptibility xtop on the 
number of colors. The instanton density (N/V) was assumed to scale as (N/V) ~ N c . The dashed 
lines show fits of the form a\N c + a2 (for ('tpiji) and N/V) and 02 + a^/Nc (for Xtop)- 
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FIG. 9: Spectrum of the Dirac operator in a pure gauge instanton ensemble for different numbers of 
colors. The eigenvalues are given in units of the QCD scale parameter. The results were obtained 
using numerical simulations with N = 128 instantons. 
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FIG. 10: Local chirality distribution function in a pure gauge instanton ensemble for different 
numbers of colors. The local chirality X is defined in the text. The results were obtained using 
numerical simulations with N = 128 instantons. 
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FIG. 11: Comparison of perturbative and instanton contributions to connected and disconnected 
correlation functions. The perturbative contribution to the disconnected correlator Fig. b is sup- 
pressed by a factor 1/N C compared to the connected correlator Fig. a. The single instanton con- 
tribution to the two correlation functions, shown in Figs, c and d, is the same, up to a sign. 
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FIG. 12: Constituent quark mass Mq, pion mass m^, and pion decay constant f n as a function 
of the number of colors N c . All quantities are given in units of the QCD scale parameter A. The 
results shown in this figure were obtained using the mean field approximation. 
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FIG. 13: Correlation functions in the pion, rho meson, and rj meson channel. The correlators 
are shown as a function of the distance in units of the inverse scale parameter. The correlation 
functions are normalized to free field behavior, Ilo(x) ~ N c /x 6 . The results shown in this figure 
were obtained using the mean field approximation. 
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FIG. 14: Meson correlation function in a pure gauge instanton ensemble for different numbers of 
colors. We show the correlation function of the pion, the rho meson, and the rj' meson normalized 
to the corresponding free correlation functions. 
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FIG. 15: Masses of the pion, the rho meson, and the rj' meson extracted from meson correlation 
function in a pure gauge instanton ensemble for different numbers of colors. The results were 
converted to physical units using A = 0.2 GeV. The quark mass was chosen to be m q = 0.2A = 40 
MeV. The dashed lines show fits of the form a± + a^j 'N c (for and m 2 p ) and a^/Nc (for m^). 
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FIG. 16: Expected behavior of the quenched pion and rj correlation functions in the large N c limit. 
The correlation functions were computed from the spectral representation equs. ( ^J49| ) using the 
phenomenological values m T = 139 MeV, = (450 MeV) 2 , E = 1.3 GeV and m = 900 MeV 
(N c = 3). We assume that \\ ~ iV c and rag ~ l/N c . 
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